Fluidization of a vertically oscillated shallow granular layer 
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Molecular dynamics simulations are used to study fluidization of a vertically vibrated, three- 
dimensional shallow granular layer. As the container acceleration is increased above g, the granular 
temperature and root mean square particle displacement increase, gradually fluidizing the layer. For 
nearly elastic particles, or low shaking frequencies, or small layer depths, the end of the fluidization 
process is marked by an abrupt increase in the granular temperature and rms particle displacement. 
The layer is then fully fluidized since macroscopic, fluid-like phenomena such as convection rolls and 
surface waves are observed. Increasing the total dissipation (by either decreasing the restitution 
coefficient or increasing the total number of particles) decreases the increase in granular tempera- 
ture and rms particle displacement at fluidization, and shifts the increase to higher accelerations. 
Increasing the frequency also decreases the magnitude of the jump, and shifts the change to lower 
accelerations. 

PACS numbers: 45.70.-n,64.70.Dv,05.70.Fh, 07.05.Tp 
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I. NOTE FROM AUTHORS 

After further investigations, we find that the results 
for the temperature and rms displacement at low T 
depend on the functional form of the velocity depen- 
dence of the restitution coefficient. Changing the de- 
pendence from e = maxfeo, 1 — (1 — eo)(v n / y/gd) 3 ^ 4 ] to 
e = max[eo, 1 — (1 — eo)(v n / y/gd) 1/>5 ] significantly reduces 
both Tki n and Ar rms at low T. The event-driven algo- 
rithm used in the present work has been tested and found 
to work well for situations where most collisions occur at 
high velocities, but this algorithm is less accurate for sit- 
uations where many collisions occur at low velocities, as 
in the present work. Hence the results we have presented 
should be checked with an model that is more accurate 
for low velocity collisions. 



II. INTRODUCTION 

Granular materials are collections of many dissipative 
particles, which can behave as either a solid or a fluid Q. 
To maintain a fluid-like state, energy must be continu- 
ally supplied to the grains since it is dissipated during 
collisions. Often this energy input is realized by shaking 
a container vertically; the non-dimensional shaking ac- 
celeration is given by T = Air 2 Af 2 / g, where g is the ac- 
celeration due to gravity, A the shaking amplitude, and 
/ the shaking frequency. With increasing T the behavior 
of the granular layer changes from solid-like to fluid-like. 
This fluidization does not occur abruptly at a particular 
T but rather is a process that develops over a range in T. 
Even for T < 1, there is relative motion of the particles 
[3, Q , and the contact forces between particles show rich 
dynamics [H . As T is increased above unity, the particles 
start to move randomly on length scales small compared 
to their diameter. This has been measured with laser 
speckle methods, where it was also found that in deep 
layers, upper layers fluidize first and then the fluidiza- 



tion proceeds downward Q with increasing T. The onset 
of macroscopically visible motion occurs at T « 2, de- 
pending on the experimental conditions [||. 

Quasi two-dimensional geometries allow the study of 
the fluidization of vertically vibrated layers using direct 
imaging. Studies of horizontal monolayers have examined 
crystallization [7} and the coexistence of different phases 
that occur during fluidization [1, [l(| HH, In a 

vertical, initially crystalline layer a distinct jump in the 
height of the center of mass of the particles was found as 
the layer fluidized [HI]. 

We present here results from fully three-dimensional 
molecular dynamics simulations of a shallow granular 
layer. The process of fluidization is accompanied by a 
smooth increase in temperature and particle mean square 
displacement that culminates in an abrupt increase for 
high coefficient of restitution, or low shaking frequency, 
or small layer depth. We identify this abrupt change with 
the end of the fluidization process. 



III. MOLECULAR DYNAMICS SIMULATION 

The event driven molecular dynamics simulation de- 
scribed in fl4| was used for this study. The number 
of frictional, rotational spheres of diameter d used was 
varied to determine the dependence of fluidization on 
the layer depth, N. The total number of particles was 
N x 248, corresponding to a volume fraction of 0.58 at 
rest. The container was a box 15d x 15d x 300d, in the 
x,y, and z directions, respectively, with the z direction 
being opposite to gravity. The bottom plate oscillated 
sinusoidally, zp = Asin(2irft). The non-dimensional fre- 
quency, /* = fy/Nd/g, was varied from 0.08 to 0.25. Un- 
less otherwise noted, periodic boundary conditions were 
imposed in the x and y directions to avoid convection 
induced by frictional sidewalls [HI, [H, Qj], EH • Surface 
waves were prevented from forming by choosing the hor- 
izontal dimensions (15d) to be smaller than the pattern 
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wavelength, which was found to be 39c? in simulations in 
a large box for a layer with depth N = 3.6, frequency 
/* = 0.15, and T = 2.3. 

The parameters characterizing the particles are the co- 
efficient of normal restitution e (the ratio of the rela- 
tive normal velocity after collision to that velocity before 
collision), the coefficient of friction \x, and the tangen- 
tial restitution /3, which gives the change in the relative 
surface velocity [14], [l9| . For this study, \i and (3 were 
fixed at n = 0.5 and (3 = 0.35, which were the val- 
ues used to reproduce surface wave patterns observed 
experimentally in oscillated layers of lead spheres [T3 |. 
The normal restitution e depended upon the relative nor- 
mal velocity, v n , of the colliding particles according to: 
e = max[e ,l - (1 - eo)(v n /y/gd) 3/4 }, as in The 
minimum value of e, eo, was varied from 0.98 to 0.65. 

Simulations were started from random initial condi- 
tions at r = 3. After 300 plate oscillation cycles, T was 
reduced by 0.02. Particle positions and velocities were 
recorded in at the phase where the plate was at its mean 
position and moving upwards. This procedure was re- 
peated with r decreasing by 0.02 each time until the par- 
ticle collisions became too frequent for the event-driven 
algorithm to be usable 20]. This limited our investiga- 
tions to r > 1.6, but we found that the accessible range 
in r included the transition to full fluidization. 



IV. RESULTS 



cycle(Fig. [ljb)) exhibits an increase at Tf, just as found 
for T kin (FigUJa)). 

We define as fluidization the process that begins when 
r < 1 and ends at Tf, the acceleration where the abrupt 
increase in Tkin and Ar rms occurs. After the increase, 
the layer is defined as vibrofiuidized. This transition is 
history independent and non-hysteretic: it occurs at the 
same T whether decreasing T from 3.00 or increasing from 
1.84, for step sizes that were as small as 0.01. 

The granular temperature and density of vertically os- 
cillated layers depend on the phase of the plate oscillation 
because, after each impact of the layer with the plate, a 
shock propagates upward through the layer [2l|. The 
phase dependence of Tkin and Ar rms was examined for 
twenty equally spaced phases throughout a cycle, and 
all but one phase were found to yield similar behavior 
(Fig. [2]) • The one phase that yielded a different depen- 
dence on r is labeled 4 in Fig. [21 at this point in the 
cycle the plate is near its minimum height and moving 
upward, which is where the shock has the greatest ef- 
fect [2lJ. However, the shock quickly travels through the 
layer, and its effects are negligible for most of the cycle. 
We therefore recorded data at the phase when the plate 
was at its mean position and moving upward. 

Just as the vertical forcing may add a dependence on 
the phase of the driving cycle, it also creates gradients in 
the z direction. The top layers are more fluidized than 
the lower layers, as Fig. [3] illustrates. While the flu- 
idization process proceeds more rapidly for the topmost 
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The horizontal granular temperature Tkin is the kinetic 
energy of the relative random motion of the particles, 



Tkin 



2gd 
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where the averages arc taken over all the particles and 
Tkin has been made nondimcnsional by dividing it by the 
gravitational potential energy of a single particle raised 
by a height equal to its diameter. The averages are taken 
only in a horizontal plane to minimize the influence of the 
vertical shock wave 21]. An example of the dependence 
of Tki n on r is shown in Fig[lja). A sharp increase in 
Tkin is evident in the range 1.94 < T < 1.955. 

We compare Tkin with the root mean square displace- 
ment, 
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where x(t) and y(t) are the x and y positions of a particle 
at time t. To minimize the impact of particles crossing 
the periodic boundary on the measurement, the aver- 
age is computed over all particles at least 5d from the 
boundary at t = 0. The rms displacement for At = 1 



FIG. 1: Comparison of different possible order parameters: 
(a) horizontal granular temperature, Tkin, and (b) horizontal 
root mean square particle displacement in one cycle, Ar rms , 
as functions of T. Simulation parameters were eo = 0.98, 
N = 3.6, and /* = 0.15. 
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layers, it does not complete until all layers have fluidized 
©. Similar behavior was found for /* = 0.25. 



B. Dependence on System Parameters 

Changing the normal restitution coefficient has a dra- 
matic effect on the granular temperature and rms dis- 
placement, as shown in Fig. |4j The transition to the 
fully fluidized state is most pronounced for our highest 
restitution value, eo = 0.98, and becomes less obvious 
as eo is decreased; for the lowest restitution value exam- 
ined, eo = 0.65, no transition is discernible. As eo is 
decreased, the fully fluidized state is reached at higher T: 
for e = 0.98, T f = 1.95; for e = 0.85, T f m 2.3. 

The dependence of Ar rms on T and eo is compared 
for three frequencies in Fig. [5j Increasing the frequency 
makes the fluidization process more gradual and de- 
creases Tf. For example, for /* = 0.08 (with e = 0.98), 
Tf = 2.0, while if /* is increased to 0.25, Tf = 1.85. 

Increasing the layer depth leads to a less pronounced 
increase in T^n and Ar rms , and the increase is shifted to 
higher T (Fig. [6]). For N — 1.8, the end of the fluidization 
process occurs at Tf = 1.48, while for N = 5.4 the change 
is smaller and Tf ri 2. Thus, as the total dissipation of 
the system increases, either by decreasing the restitution 
coefficient or by increasing the number of collisions (by 
adding more particles), the change from the fluidizing 
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FIG. 2: Phase dependence of the (a) granular temperature 
and (b) root mean square particle displacement in one cycle 
as a function of T. The inset in (b) shows the height of the 
plate (normalized by the oscillation amplitude), zp/A, as a 
function of the phase tf of the driving cycle. The temperature 
curves for different phases in the cycle are similar except for 
phase four (open triangles), which is influenced by the shock 
wave that passes through the layer after the particles hit the 
plate [2lj . Filled circles indicate the phase used for the data 
presented in other figures. Simulation parameters were eo = 
0.98, N = 3.6, and /* = 0.15. 



TABLE I: T onset , the onset of fluid- like phenomena, is found 
to be greater than or equal to Tf, the beginning of the fully 
fluidized state. 
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to the fluidized state becomes less pronounced and Tf 
increases. 

After the layer has completed the fluidization process, 
fluid-like phenomena appear in addition to the increases 
in Tki n and Ar rms . For example, as is well known, 
surface waves appear at T 2.5, [3 [H, [23j], where 
T f fa 2 [(| . We have done some simulations in a container 
large enough to accommodate waves and have found that 
waves emerge in a thin layer at about Tf, while for a 
deeper layer, pattern onset occurs for r > Tf (see Ta- 
ble Hj). Additional simulations were made for a container 
with solid frictional lateral walls, and a single convection 
roll was found to develop at T = 2.24, the same T at 
which the fluidization process completes (Table [T|. 



C. Comparison with Experiment 

Comparing these simulations to experiment is prob- 
lematic. Measurements of the granular temperature or 
rms displacement in a three-dimensional experiment are 
challenging since the trajectories of the grains in the bulk 
of the layer can not be recorded with just a simple video 
camera. Recently speckle visibility has been developed 




FIG. 3: Height dependence of (a) T^ in and (b) Ar rms as 
functions of T, measured above the surface of the plate. Flu- 
idization completes at the same value of T for each height, 
but particles higher in the container are more mobile at lower 
P. Here e = 0.98, N = 3.6, and /* = 0.15. 



4 



to measure the granular temperature in a 3D sample [24] 
and should be utilized to fully explore this process. 

Umbanhowar and Swinney investigated the transition 
to surface waves for systems with similar layer depths 
and frequency using the pressure exerted by the layer 
on the bottom plate. They therefore recorded data at 
a different point in the driving cycle than was used in 
this work. In addition, the fluidization process may have 
completed at a lower T value than they report for the 
pattern transition as in our own work (Table [I]) [22j ■ 

Two groups have studied fluidization in deep layers 
(N = 15d in [|| and N > 20d in [5]) and have reported 
that the fluidization of the layer begins at the top and 
proceeds downward as T is increased. Evidence of this 
top down fluidization can be seen in Fig. [3] Particles 
higher in the container are in a more fluidized, but flu- 
idization is completed at the same value of T for particles 
at all heights studied. 

V. DISCUSSION 

We have shown that in shallow layers of nearly elastic 
particles, abrupt increases of both the granular temper- 
ature and the root mean square displacement indicate 
that the process of fluidization is completed. If the dissi- 



pation in the system is increased by adding more particles 
or decreasing the coefficient of restitution, the increases 
in granular temperature and rms displacement at full flu- 
idization become more and more gradual. 

While the granular temperature and root mean square 
displacement are good indicators of the onset of the vi- 
brofluidized state, the question remains if there exists a 
single control parameter for the fluidization of a verti- 
cally oscillated granular sample. The non-dimensional 
acceleration T is not satisfactory since layers with differ- 
ent characteristics fluidize at different values of T. The 
value of the temperature Tki n depends on the phase of 
the cycle, and both the rms displacement and Tki n at 
the onset of the vibrofluidized state depend on granu- 
lar parameters such as coefficient of restitution and layer 
depth. Thus they do not make good measures. Others 
have suggested vq = A{2-kJ), but this measure does not 
collapse our data since we find a different dependence on 
/* than that reported in [i~3| |. We hope that our results 
will inspire further work on the process of fluidization in 
three-dimensional granular layers. 

The authors thank William McCormick and Chris 
Kriille for helpful discussions. This work is supported 
by the Robert A. Welch Foundation grant F-0805. 
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FIG. 4: (a) Linear and (b) logarithmic plots of the granu- 
lar temperature, and (c) root mean square displacement as 
a function of Y for different values of the normal restitution 
coefficient. The fluidization of the layer is more abrupt for 
particles that are nearly elastic. Simulation parameters were 
/* = 0.15 and N = 3.6. 
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FIG. 5: The effect of varying the frequency on the F depen- 
dence of Ar rms . The change signaling the end of the fluidiza- 
tion process becomes more pronounced as the frequency is 
decreased. Here N = 3.6. 
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FIG. 6: The effect of varying the layer depth N on the 
r dependence of Tkin and Ar rms . As the layer becomes 
thicker, the change signaling the end of the fluidization pro- 
cess becomes less pronounced. Simulation parameters were 
/* = 0.25 and e = 0.98. 
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